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ABSTRACT 



Aircraft m fin ■hnum fTy.ing; speed:, as determined:by actual flight test, 
is published in aircraft handbooks for pilot’ guidance; The test flight 
results are used to determine and confirm takeroff-and rlanding speeds, 
field TengthSr left-hand: portion of the maneuvering envelopes (V-n 
diagram);, etc. Beterminatd'an oT the absoTuteeminimum flying speed of 
an aircraft on the other* hand:,, has: not been of* prime. importance in 
flight test. 

In the present anaTysn’s: digital simulation allowed the systematic 
study of not only the minimum flying speed, assdefined .by Federal Aviation 
Regulations but also the absoTute minimum flyingrspeed .attainable in 
steady, unaccelerated flight*. The. study included:such effects as 
deceleration rate,, rate orf change of elevator- anglei .aircraft weight 
and pitch moment of inertte.. 

It was found for an ctssumed Tight-weight'.fighter aircraft that the 
absolute minimum flying’ spee± was: approximately 20 knots less than the 
FAR minimum flying speed.. Moreover the FAR- minimum flying speeds 
tended to be quite sensitive to rate of change of elevator angle. 
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E.. IWTRDDUCTION 



Kimimum flying speed and or- sdalTing: speed: is. defined by, Federal 
Aviation Regulations (FAR) [iTef'.. IT as: ttieinrinimunr steady speed rat • 
which the aircraft remains crmtral Table: . The: Naval TestrPi lot .School 
Flight Test Manual [Ref.. Zl defiries: minimum" flyings speed, as :the . minimum 
steady airspeed attainable in unaccHlerated-flight:or''the:minimura 
usable airspeed. The actual determination of' mini mum flyi ng .speed :is 
always accomplished through fTrght. tests: . Methods: for 'theoretical ' 
predictions are empirfcal,. bashed, upon V'/indrtunneT results, .experience, 
etc. 

FAR specified that stalls^ would be^ demonstrated: by. trimming :the 
aircraft at 130 percent of the estimated minimum flying, speed .and , 
decelerating at a constant rate- until’ a_mihiinum speed. was lobtained. 

The certified stall speed would correspond. to a: deceleration rate of 
1 knot/sec. To facilitate the test, pilot's: task, .Ref; .1 and .ZIdefined 
aircraft characteristics and parameters^ indicative:of 'minimum flying 
speed. Typical indications were:- longitudinal , directional or lateral 
divergence, excessive altitude loss, loss; of control /effectiveness , etc. 
Data analysis methods and restrictions' on demonstration technique, 
other than rate of change of airspeed,, were: not: specified by Ref. 1. 

Additional flight test techniques have; been proposed and .evaluated 
by the National Aeronautics and Space Administration [RefsVS.’and 4]. 
These studies, through flight tests:,, evaluated" three; demonstration 
techniques: FAR, T-g break and constant-raterof-climb. Three methods 

of data-analysis were used with the FAR’ demonstration technique .to 



determine minimum fTying speed:; cmnsdant: deceil erati on-, ; and:C|l methods. 

The minimum flying speeds' obtafned by theE thre^ techniques: wereicorrected 
to a deceleration rate of T. knot/ sec id aonardance: witb Ref: . T. . 

The study described herein appTied. numericaT methods: to: simulatec 
flight test data in an attempt to define an absolute: minimum* flying: 
speed. Five data-analysis methods were used.. The speed. was: an- absolute . 
minimum obtainable for steady,. unacceTerated flight'. when- the: FAR- 
correction to a 1 knot/sec.. deceTeration' rate was: disallowed: . The: 
absolute minimum ■Flying speeds were, compared with' FAR- certified:minimum 
flying speeds, i.e., corrected to T knot/sec.. deceleration- rate; . 

“Flight test data" for the- F-34A,, a smalT snngle^engine: jettffghter, . 
v/as obtained from a computer program of the non-l ihear equationsrof 
aircraft motion. FORTRAN IV language and the: Naval' Postgraduate :School ‘ 
IBM 360 digital computer were utiTrzed.. The:- computer^ program provided, 
numerous options For trimming and flying the. aircraft*,, ihcluding. 
variations of thrust, weight,, rate of change of* elevator angle. and: 
flight orientation, i.e.., cTimbfng,. level' or- descending- fli ght: 

For this study, the aircraft was trimmed in* steady, .descending ; 
flight at 130 percent of the estimated minimum flying speedand decelera- 
ted by decreasing the elevator angle at a constant rate. Two values, 
for thrust were used to trim and stall the aircraft; T/W s.0^02.and 
T/W = 0.11. 

The effects on minimum flying speed of thrust, aircraft.weight, 
rate of change of elevator angle and aircraft moment: of' inertia :about 
the lateral axis were also studied. 
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Hie purpose of tJifs study, was:- toe invest! gate: the; relationships of 
the aircraft mlnTmum fTyitigf speeds: as; determined: by. the;different data- 



analysfs methods.. 



O.. EXPERIMENT7\L methods: 



A. BASIC APPROACH 

The study of minimum fTying s:peeit wets; restritted: to; aircraft'.motion 
in the XZ plane: i.e.* to aircraft: motion ailbng: the' X' and: Zlaxesiand. 

pitching moments about the Y axrs:-.. I!t wass assumed: that: no: adverse: 
lateral or directional traits occiirre± while: dece3era±ing:to:the: 
minimum flying speed. 

Non-linear equations of aircraft: motibrt with’' threec-degreesfof-'- 
freedom v/ere derived fron EuTer' s equations: of'" motion,, Appendix; A: . 
Extensive use was made af Etki'n 5c]'.. The;- equations. wereeset:up. 

for input of non-linear aerodynamic data,, by means: of* a: table; look-up. 
procedure. 

A computer program was wri tten irr FORTRAN IV. ranguage:fdr’'the;Naval ' 

Postgraduate School IBM 360 digital' computer,. Appendix. Bl. Theiprogram 

employed a fourth-order Runge-Kutta algorithm in' sol ving the.equations . 

of motion. Real-time was used in eval uating: the; aircraft:motion''from 

approximately 1.3 to Typical' time: histories, of* the 

aircraft parameters were plotted in Figure T. Rgure. 2-contains'Summary 

plots of V as a function of deceleration rate for different: rates 
min 

of change of elevator angle. The thrust to weight ration (T/W) was:0.02. 

B. TEST AIRCRAFT 

The F-94A, a small single-engine jet fighter ,. was: uti 1 i zed as the 
test aircraft. 

Reference aerodynamic data were obtained, from Bl'akelock [Ref: 5] 'in 
the form of linear stability derivatives for alT terms except'those 






Figure 1. Typical Time Histories of Aircraft Parameters 
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Figure 1. (Cont'd) 
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depending upon angle of attack-.. Thee vea'iation of- theecotnpleteeaircraft 

C^, Cq and Cj(| was estimated: from the available data: for aecenter of 

gravity located at the Q.ZSZJT (wirg-mearr aerodynamic:cbord) station 

and an airplane in the landing cmnfiguration. Although theevariations 

of C^, Cq and Cj^ with angle of a-ttack,, a,, were estimated for zero 

elevator deflection* the use of the linear- stability derivatives such 
C C 

as L. and M- allowed trinming thee aircraft to any initial flight 

0_ Q.. 

e e 

conditions. The C,, vs., a curve for 6' = 0' was selectedcfdr a;trim 

C. of about I AS with the sdalT break being very gradual ; .which is 
*'max 

indicative of a progressive fTov/ separation starting- at:the iwing trail ing 
edge. 

Aircraft reference data cmd physdtaT characteristicsiareilisted in 
Table I. 

The choice of afrcra-ft fTfght traits: are arbitrary, and :may, seem 
restrictive to a person who is solely concerned about:a:particular 
aircraft by virtue of personal experiences. However, . i t'.sbould be 
recognized that the methods and qualities, described- in- this. analysis 
are typical and may be modified tc a: particular application by. changing 
the input data. 

C. FLIGHT PROCEDURES 

All "test flights" were conducted with the F-94A configurated for a 
landing approach. The aircraft was trimmed for steady descending flight 
at 130 percent of the estimated minimum flying speed. Aircraft thrust, 
vreight and moment of inertia about the lateral axis. were. the .only 
physical characteristics changed during- the "test flights". 
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TABLE I 



oE. PhystcaT Data 



Parameter 


Value 




239'ft^.- 




6£4-:ft: 



12;359:ib 



V 

((Tinear region) 
a. 


26i'545:slug-ft 

5£2Z>ad’^ 


\ 

e 


0C433rad“^ 


Cj^. (n.itiear region) 

CL 


-0:406 rad’ ^ 


\ 


-0:88:rad’^ 


si. 
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-4r268 rad"^ 




-8.165 rad'^ 


c 


25:2% c 
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MBLE E (Carrt'd.) 



b. Tabular Data (complete aircraft, 6^ = 0^ 



(rad) 






"'o: 


0.10472 


0»55Z 


-(LDUZr 


o: 09405 


0.13963 


G»738 


-(Lono: 


o:ii5o. 


0.17453 


a.9ZT 


-<L0335; 


o:i4oo: 


0.20944 


Ti»Q93- 


-0104455 


o: 1695 5 


0.24435 


T:..232 


-OC. 05855 


o:203o: 


0.27925 


T..346: 


-05.0740: 


0C2490: 


0.31416 


T..430 


-0109.10: 


o:3i2o: 


0.34907 


T>484 


-oiiizo: 


0'.3640: 


0.38397 


T»5T3 


-0C..132O: 


0C40505 


0.41888 


T.520 


-0C.152D: 


014550 5 


0.45379 


T.500 


-01172D‘ 


014930' 


0.48869 


T>46Z 


-OC. 20755 


0153505 


0.52360 


1..41.5 


-a. 2400: 


0157205 


0.55851 


T.350 


-0C.Z795: 


016100' 
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Deceleration from the trrnr cxmxii'ticrrr was: accompl'i shed: by.. decreasing ; 
the elevator angle, to stmlate aft cnntroil' stict moyement: . Available -: 
time histories, [Ref. frrdrrcated: elevator- angle: to:be:a:neac 'linear 
function of real-time during- tJie deceleration maneuver;. Therefore^ a: 
constant rate of change of elevator ang.lewas: chosen- for this .'Study,. 

Rate of change of elevator ang^le was: varied: between--Gj05.and:-K15.' 
deg/sec. for the different maneuvers:.. 

The effects of thrust were studfed. by making, two: "test:flights“ for 
each elevator angle schedule.. Two thrust values: were: used: fdr'-the . 
“flights" which gave thrust to wefght ratios: of- 01021 and. Oil K . 

Five "flights" were made at -0;.28ff deg/ sec:, rate: of-' change :of-' 
elevator angle while varying aircraft wei ght to: study thereffects:of ' 
aircraft weight on minimum flying speed.. The aircraft weight'.was : 
varied from 12,355 pounds to Iffj.JSS poundsMh- incrementsiof ''1000 Ipounds . 

“Flights" were made at representative rates' of- change* of -'elevator 
angle to study the effects erf aircraft moment ot inertia about: the: 
lateral axis, on minfmum flying speed., Values: 251 percent' above ..and : 
below the true aircraft moment of inertia were: evaluated. 

Data from the original "test: flight" indicated the possible-:exist- 
ance of an optimum elevator schedule in determining, minimum flying 
speed. An exponential elevator schedule was studied and. compared to 
the linear elevator schedule.. The aircraft v/as decelerated using a 
linear schedule until the deceleration ra-te" fell below a. preselected 
rate. At this point the elevator schedule. waS; changed to an exponential 
schedule until minimum flying speed was obtained. 
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D. DATA-ANALYSIS METHODS 



Five data-anaTysfs methads were:; used tix dedermi ne-; mi nmum fTyingj 
speed from the tv;o mafn sets af "fT.ight: tesd!' dates . Theseemethods.were 

1 . Constant DeceTerattorr Rate 

The linear elevator angle: S3cheduTe resulted: in- near constant: 

deceleration rates by the test afrcraft.. Mirtimum flying:speed:was : 

defined as the actual nrrnimum speed obtaitied during, the deceleration. 

For consistency, a deceleratioir rate was defined' as.- tbe:s3ope:of''a: 

straight line drawn from V to T..T, V This: deceleration rate:was ' 
^ mtn min 

used in both the Cj^ and methods.. Plots: were made-: of-' minimum flying 
speed as a function of deceTeratiou rate and: minimunr flying-speed _as .a 
function of rate of change of elevator- angle: for- the: two. values : of ■' 
thrust. 

The constant decelerattoiT rate method was: used: to: analyze . data ; 
when the effects of aircraft weight and moment, of inertia: on-minimum 
flying speed were studied. 

2. Methods 

was defined as the aircraft left coefficient: independent: 
of normal flight path acceleration: 

W 



The assumption of aircraft weight and aircraft lift being equivalent 

throughout the maneuver corresponds to considering C|^ as. being, an 

apparent lift coefficient rather than a true value. Minimum flying 

speed was defined by the Cj^ method as the speed of the. aircraft when 

C|* was at a maximum value.. Plots were madeo.f Ci* as: a: function of 

^max 

21 ' 



deceleratton rate,. ^'C' ace e- furrctibn' of' deceleration •• rate cand :'^C ' . 

maX’ mi n 

as a functTQiT of rate of Ghange- Qf-'eieyator angle. 

3. C|^ Method 

C|^ was defined: as beiiig: the true: aircraft' 1 1ft '.coefficient which 
compensated for effects of normal' flight*. path acceleration: 






W,.a: . 



qs: 



K'in.Tmum 'fTytng speed was defined by, thecC^ method asctherspeed .of the 

C ' 

a 'ire raft when C, was at a raa'xiinum value. . Plots, were imade :of L' as 
L max 

a function of deceleration' rate,. as: a function- of ''deceleration 

mm 

rate and as a. function of rate of change of*elevator angle, 

mm ^ ^ 

4. Altitude Break Method 

The aircraft aitttude trace had' a relatively, 1 1 hear character- 
istics slope during the deceleration maneuver. With- theconset of stall 
the characteristic slope would increase significantly in- value.' The 
minimum flying speed was de-fined as. the speed of the. aircraft where the 
linear characteri Stic slope of the altitude trace could .no longer be 
maintained. The characteristic slope increase was. indicative of 
increased altitude loss. 

A deceleration rate was defined as the slope of a straight line 

drawn from '^h„. to 1 ..1 '^h . . Plots were made of as a function 

mm mm mm 

of rate of change of elevator angle.. 

5. 1-g Break Method 

The minimum flying speed defined by this method was the last 
speed at which the aircraft could: maintain 1-g flight;. This method 
proved unreliable in data reduction for an aircraft trimmed in steady, 
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descending flight and was thus not pursued further-.. For- ttiiss method', 
an aircraft must be decelerated from an mi tTatT cTiinbr. condition. 
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nn£.. RESULTS ANDL DISCUSSION 



A. ABSOLUTE MINIMUM FLYING SPEED 

1. Constant DeceTerattan Rate Method- 

The minfinuni fTytng speeds de.termined;were=pI6tted.in Figure 3 
as a function of deceTeratiorr for thrus1:-to'-weight*-ratiosp0.'02 and 0.11. 
For a thrust-to-wefght ratio aP D.OZ the. absoJuteitninimura flying speed 
was S3.7 knots and occurred at c6 deceleration rate of^2160Iknots/sec. 

A thrust-to-wefght ratio of Q..1T decreased, the absoluteTminimura flying 
speed to 80.6 knots with na appreciabrTe change, in' deceleration-rate. 
Figure 2 also indicated a b salute minimum flying; speed:;v;as:sensitive to 
variation in optimum deceleration rate;.. The; addition of 'thrust 'reduced 
this sensitivity somewhat.. 

Figure 4 indicated absoTute minimum- flying: speed .was .'relatively 
insensitive to rate of change of elevator.. Ra.tes'between- -0.7 "and -1.1 
deg/sec. were able to give a close approximation' of -'absolute minimum 
flying speed. 

• 2. Method 

Maximum values of were plotted in Figure. 5: as a function of 

deceleration rate. . was plotted as a function of deceleration rate 

mtn 

in Figure 6 and as a function of rate of change of-'elevator angle in 
Figure 7. A C|^ of 1..815 was the maximum attainable in steady, 
unaccelerated flight.. The absolute minimum flying speed was 83.5 'knots. 
The thrust-to-v/eight ratio of ff.lT decreased: the absolute minimum flying 
speed three knots with minimaT effects; on deceleration'rate. . was 

fairly sensitive to deceleration rate and. insensitive to.rate of change 
of elevator angle. 
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Figure 3. Minimum Flying Speed as a Function of 
Deceleration Rate. 
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Figure 4. Minimum Flying Speed as* a' Function of- Rate of 
Change of Elevator Angle. 
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Figure 5. Apparent Maximum Lift Coefficient as a Function 
of Deceleration Rate. 
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Figure 6. Minimum Hying Speed as a Function of 
Deceleration Rate. 
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Figure 7. Minimum Flying Speed as a Function of Rate 
of Change of Elevator Angle. 
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3. Method 



Absolute minimum flying speed was 88 knots: for' ee thrustrto:— 
weight ratio of 0.02, Figures 8 and 9.. Irrcreasing. thrusii decreased: 
this speed 2.5 knots. was sensitive ta hath deceleration' rate: and: 
rate of change of elevator angle. Figure TQ indicated; the abso:rute: 
maximum C|^ value v/as insensitive to deceTeratfon rate and increased: 
thrust. 

4. Altitude Break Method 

V 

*^min plotted as a function of deceleration* rate: in- Figure: 

11 and as a function of rate of change of elevator- angle: in Figure: 121 . 

The absolute minimum flying speed for a thrust-ta-v-;eig,ht ratio: of01021 

was 83.5 knots at a deceleration rate of -£..6 knots/ sec.. Rgure; IV 

indicated that two values for existed for* one: val ue: of deceleration 

mm 

rate between -2.6 and -3.76 kts/sec. Figure T£ mdicated wass 
insensitive to rate of change of elevator angle.. 

The altitude break method was dependent on graphical' interpretaf- 
tion for data. This accounts for the smaTT amount of data scatter* in-^ 
Figures 11 and 12. 

5. Comparison of Methods 

The general relationships of absolute mini-mum flying speed as. 
determined by the various methods are indicated in Figure 13.. 

The constant deceleration, Cj^ and altitude break methods defined 
the absolute minimum flying speed to be 83 knots, Rgure 13. The. 
method absolute minimum flying speed was 83 kno-ts. Increased, thrust* 
decreased the absolute minimum flying speed defined* by the four methods, 
about three knots. 
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Figure 8. Minimum Flying Speed as a Function of Deceleration 
Rate. 
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Figure 9. Minimum Flytng^ Speed as a- Function of Rate of 
Change of Elevator Angle. 
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Figure 10. True Maximum Lift Coefficient as a Function 
of Deceleration Rate. 
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Figure 11. Minimum Flying Speed as a Function of Deceleration 
Rate. 
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Figure 12. Minimum Flying Speed as a Function of Rate 
of Change of Elevator Angle. 
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Figure 13. Comparison of Aircraft Parameters at Absolute 
Minimum Flying Speed T/W < 0.02. 
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The deceleratfon rate was apprcixiTii2Et^y-216c knots/see;at:thee 
absolute miniraura flying speed far aTT four methods:. . S3gnificant*.changes 
in deceleration rate due ta increased thrust wass observed: only, in'the; 
altitude break method. 

Minimum flying speed determined by the: method: occurred :at'. 

the smallest rate of change of eTevatar angle,. -014^ degs/sec; . Minimum 
flying speed determined by the constant deceleration,; and:altltude; 
break methods occurred at -Q..9 deg/ sec., rate of-' changer ofelevator- angle 
Absolute minimum flying speed defined by the C‘|^ method, occurred .at :the -: 
smallest elevator angle, approximately -9 degrees:,; v/hereas: approximately 
-13 degrees elevator angle was required for the other methods: . 

The angle of attack at the absolute minimum flyi ng' speed :was ~ 
approximately 28 degrees for the constant deceTerati on , . and.altitude 
break methods, and considerably less at 2T1‘ degrees: for the: method. . 
Thrust effects on angle of attack at the ahso.Tute: minimum flying-speed; 
were insignificant (less than 0.5 degrees) for alT methods. . 

Elevator angles, lift coefficients and angle; of attacks, at: 
minimum flying speed/stall speed were less than,, in some: cases; .the: 
maximum static values defined for the aircraft. Possible: explanations . 
are: 

1. Rate of change of elevator angle was. programmed: to become 

zero two seconds after was obtained.. For small . values: of elevator 

mm 

rate, occurred much sooner than the maximum obtained values.of -C^, 

a and 6^. Therefore, elevator angle was held constant until the-:end of 

the maneuver, restricting the aircraft from obtaining the maximum values 

of C, , a and 6„. 

L e 






2. The variatfons of drag cxiefftcient: w11:h angiee of ''attack- used: 
in the table look up were represerrtattve vaTues and not: necessarily true: 
values for the F-94A. Reduced drag coefftcients would: therefore: effect: 
the occurrence of minimum flying speed. 

3. In defense of the low values of C a and 6' it:must:be-: 

L e: 

remembered that this is a dynamic study of stalT tran'-tsi. Actual' flight: 
test data, [Ref. 4], indicates angle of attack is: less: than" static: stall 
angle of attack. 

B. FAR MINIMUM FLYING SPEED 

The minimum flying speed defined by Federal' Aviation Regulations: 
corresponds to the minimum speed obtained at a: one knot/sec: . decel era-- 
tion rate. The FAR minimum flying speed determined by the. constant: 
deceleration, and C|^ methods was approximately 105' knots. The: 
altitude break method's FAR minfmum flying; speed, was: 109.- knots; . Increased 
thrust decreased the altitude break minimum speed two: knots: and decreased 
the minimum speed defined by the remaining me'thods' Less: than one: knot; 

The rate of change of elevator was - O...OS deg/sec., at: minimum flying, 
speed determined by the constant deceleration rate,. C'|^ and: methods. 

The altitude break method rate of change of elevator angle. was -0.05 
deg/sec. 

Angle of attack was approximately 15 degrees and elevator angle -4.8 
degrees for all methods. 

C. COMPARISON OF MINIMUM FLYING SPEEDS 

The absolute minimum flying speeds were categorically 2D knots: less 
than those defined by FAR one knot/sec. deceleration rate,. FigureT4. 
Thrust had more effect on the absolute minimum flying speed whereas rate 
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Figure 14. Comparison of Aircraft Parameters at FAR Minimum 
Flying Speed. 
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of change of elevator angle had mare effect an the FAR mitiitnum flying: 
speeds. Angle of attack and elevator angle were sfgniffcantJy greater 
with the absolute minimum flying speeds. 

D. SENSITIVITY STUDIES 

1 . Aircraft Weight 

For a representative rate of change of elevator- angle,, Rgure: 

15 indicates minimum flying speed was a linear function at aircraft,, 
weight. 

From the definition: 

W . 

hnax g^S 

minimum flying speed is seen to be a function of the one-half power- of 
weight. Actual flight test data indicates minimum flying speed: is: ai 
linear function of aircraft weight [Ref. 4].. 

2. Aircraft Pitch Moment of Inertia 

When the aircraft pitch moment of inertia wac increased or- 
decreased 25 percent, there were no significant effects on minimunr 
flying speeds or other aircraft parameters. This conclusicn was. not' 
based on a systematic study but rather a brief look at three data 
runs. The trends were present but a complete study is warranted. 

3. Exponential Elevator Rate 

Comparative plots of linear and linear plus exponential e.levator 
rates were made in Figure 16 for initial linear rates of -0.06, and 
-0.115 deg/sec. The exponential rate was e^ where x was increased 0.04 
each second. obtained by the linear plus exponential elevator rate: 

was significantly less than that obtained by the linear elevator rate. 
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Figure 15. Minimum Flying Speed as a Function of Aircraft Weight. 
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Plus Exponential Rate pf change qf Elevatqj" Angle. 



TBiese results;,, though corTy. et samp! ing of- data: , . tend :to support the. 
existence of an optiinunr el'evcEtor' schedule in determining :mini mum flying 
speed. 
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lY. CQNCXUSrOIlS ATIE FECDMMENDAnONS: 



Digital simuTatian has great patemticel' iir the: eyal'uation'of-minimum 
flying speed. The options avatTah-Tc are practically, unlimited; 

The determination of ahsoTute mitiimum flying; speed; and: aircraft', 
behavior with digital simuTattarr wiTV enatiTe the. engineer to: better 
predict the aircraft charactertstfcs: encmintered' during: test'.fl ight: 

The test pilot's job will be si'mpTffied. somewhat: if' he; has 'knowledge :of 
general aircraft characterfsttce prior toe fl ight.. 

The operational pilot's performance should be: enhanced: by having 
an insight into his aircraft's performance beyond; the: certified stall ' 
speed. 

The computer program coded for this study may be; used for many, 
fields of study as indicated by Ref".. 7,. where it: was; used to:make;a: 
stability analysis in the vicinity erf minimum flying speed; 
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RON-LINEAR EQUATIDUS' OF AIRCRAPr MOTION' 

I. Euler's Equations of Aircraft Mbti on. 

The aircraft was considered. aS: ee ri^'dibody. with aiCartesian 
coordinate system fixed at the center' of mass:,., FigureeUr. The ; frame of 
reference Cxyz, a body axis system,, moved with' the: airccaft:and :Cx was 
fixed to the longitudinal axis of the aircraft.. The:-resultant external 
forces and the moments of these forces' about: the. aircraff:center of mass, 
when referred to the frame of reference: Cxyzq,, where: thei-vector equations 
of motion: _ 

_ _ 

F = nr + m oi x 

The scalar components of these equations are: Eul er 's: equations .of ■' 
aircraft motion. 

II. Aircraft Orientation and Axes. 

Reference frame Cxyz could not be used, to describe the position and 
orientation of the aircraft as it was fixed, to the center of mass -and 
moved with the aircraft. An earth fixed: reference frame, assumming 
the earth's rotation was negligible,, was- chosen to describe the aircraft 
motion. 

Euler's equations of motion are valid for any orthogonal ‘ reference 
frame fixed to the aircraft center of mass... To; simpl ify. the .equations 
of motion and expressions for aerodynamic: forces,. the: reference .frame. 



Cx'y'z', was chosen as stahtTity. aExes,, i -es..,. Cx-'' [rodhted: in' tbeedireetlon-' 
of motion of the aircraft center of mass:,, Fi'gure: IT.. 

III. Assumptions 

In addition to the assumptions previmisTy'marie:-defining;airccaff: 
orientation and axes it was assumed that aircraft. mot-ion- v/as: restricted: 
to the Cxy, longitudinal plane.. This assumptioir equated; the: following: 
quantities to zero: L, N, P,. R„ t,. V. and Y-.. Cxy. was; also: assumed 'to- 
be a plane of symmetry which equated the products; of-" inertia, , DIand-Fj . 
to zero . 



IV. Mon-Dimensional ization 

The external aerodynamic farces arid moments; were: expressed. as: 
functions of dimensionless force and moment caefficients.,. dynamiccpressure, 
aircraft planform wing area and a reference Vength.. Theseiequations; . 
in scalar form, were: 



X 

Z 



" ''c s: 

o ''I s: 



When these values were substituted into the scalar resultant external 
force and moment equations the following equations were obtained: 

P ''c ^ 

U = ^ ^ - g[ sin a - QW 

2m 

p s c 

W = + g COS' 0 + QU 

2m 
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BODY AXES 



STABILITY AXES 




Figure T7. Ai rcraf t Axes . 
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V. Aerodynamic Force and Moment Coefficients' 
From Figure 18, the force coeffici ents: v;ere: 



= Gj- + Cj^ ^ITT a. - C^- CDS; a. 



= -Cj^ cos a - Cg sin OL 



The pitching moment and Tfft coefficients: consisted-of*! inear portions. 
The non-linearity was present in C^(a) and C|j^(a) ,.which were;non-1inear 
curves of Cj^ and versus .. Complete expressions: for 'thectotal ’ 
pitching moment and Tfft coefficients were: 



VI. Final Equations 

Angle of attack appeared in the pitching moment equation as a 
derivative. To obtain a set of simultaneous differential equations that 
could be solved by one of the methods available, . another equation was 
added defining the rate of change of angle of attack. The.three-degrees- 
of -freedom, non-linear differential equations: of aircraft motion programmed 
for the Naval Postgraduate School' IBM 360 dlgitaf computer were: 
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Figure 18. Aerodynamic Force and Moment Coefficients. 
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APPENDIX B: 



DEa/aOPMENT OF- COMPUTEB- PROGRAf^ 

I. COMPUTER FACILITY 

liliie arf'gtnal iti tent, was: to: program theeequatlons :of 'motion on the 
AerorrauttcaT Engineering: OepartmentTEAL 580 1 ana log computer. Problems 
wHth, insufficient equipment toe progranr theeequations ;of 'motion on the 
analog computer necessitated the:- use of-'digital 'simulation. 

The Naval Postgraduate School'' s: IBM 360;digital computer offered 
many advantages and was: chosen- as: the; research tooK , Several, major 
advantages were:; the Targe amount: of 'corezavailable, .dimensioned 
qqantities could he used,, plot subroutinesiwere :available :for time 
fcfistory studies and. numerous methods wereeavailable .for 'solving 
simuTitaneous non-Tihear- differential' equations on a digital computer. 

II. COMPUTER PROGRAM 

®f the many one step' procedures: that'.were.available cfo'r solving 
systems of simuTtaneous flrs.trorder- ordinary, differential equations, 
a fourth-order Runge-hutta algorithm was; chosen, [Ref;‘ 8]j This method 
offered simplicity,, accuracy and; short computer time. The Riinge-Kutta 
procedures were programmed into subroutine RUNKUTI RUNKUT 'called 
function subroutines FI, F2, F3, F4 and F5 to calculate the derivative 
terms at time T(I) ,. T(I')+HI/Z and.T(I)-HHI. The derivative terms were 
multiplied by the appropriate. weight: factors and. added to the variables 
calculated at time 1(1) which gave. values for the. variables at time 
TCI+I). 
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Function, subroutine FT,. caETTert saibroutinee SPLItl toccalcalateethe : 
non-linearities in C^,, amrct Cjj.. SPUN wass obtained: from' thecsctentific 
subroutine package of the Navel Pdcsrtgraduate:- School* IBM’ SSOIdigital* 
computer and modified for this program.. Non-T *1 near tables of 'C|^, and 
Cq as a function of angle of attack v/ere: read: into: SPUN'. . When -supplied 
values of angle of attack,. SPLEN. used a cubic: spl ine: function' to provide 
interpolated values of C^^,. amd Cj^.. 

Initial aircraft trim conditions: U.,. W,. 0-,. S’,, a. and:6^.were:supplied 
to RUNKUT at time zero for initi'a-tidn of the: step, procedure. . These 
values were calculated in Subroutine TRIM' from static:fdrce:and .moment 
equations. 

Control of the aircraft was accomplished: through changes?in elevator 
angle and thrust which simulated control* stick and. throttle ;movements. 

The options available were seTected: through' the: integersion’- 
the first data card. 

III. TRUNCATION ERROR AND STEF-SIZE ANALYSIS: 

Error analysis by rigorous mathematical* derivationssarecvirtually 
impossible to implement for higher-order Runge-Kutta algorithms .for 
systems of differential equations. The truncation error: 

^t " If ^^n+1,2 "'-^n+1,1^ 

for a first-order ordinary differential equation was. applied to the 
system of differential equations as a step-size control mechanism, [Ref. 8]. 
Step-size and truncation errors are l isted in Table. If 
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TaliTe El 



Local Truncatfan Errors 



ep-Size 


Error 


0.2 


- 2 ?. 32 : X lo:"^ 


0.1 


-3?. ST X IIX’^' 
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0.025 
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Tl-E INITIAL CONDITIONS ARE READ INTO THE PROGRAM 
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